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Abstract 

Spatially homogeneous field theories are studied in the framework of dy- 
namical system theory. In particular we consider a model of inflationary cos- 
mology and a Yang-Mills-Higgs system. We discuss also the role of quantum 
chaos and its application to field theories. 
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I Introduction 

Quantum field theory offers a wide variety of applications, in particular for 
condensed matter^ and elementary particle physics^. Field theoretic ideas 
also reach for the cosmos through the development of the inflationary scenario 
- a speculative, but completely physical analysis of the early universe, which 
appears to be consistent with available observations^. 

In the last years there has been much interest in the chaotic behaviour of 
field theories^"^. In this paper we discuss and extend our recent results^^"^''' 
on instability and chaos in classical and quantum field theory. In Section 2 
we show how spatially homogeneous field theories can be studied by using 
the dynamical system theory and we introduce some basic definitions for the 
regular and chaotic dynamics of classical and quantum systems. In Section 3 
we analyze the local stability of a inflationary scalar field minimally coupled 
to gravity and its point attractors in the phase space. The value of the scalar 
field in the vacuum is a bifurcation parameter and we discuss the existence of 
a stable limit cycle. Finally, in Section 4 we study the spatially homogenous 
SU(2) Yang-Mills-Higgs system. We show that for this system a classical 
order-chaos transition occurs both in classical and quantum mechanics. 

II Field theories as dynamical systems 

In this Section we introduce some basic ideas of the dynamical system theory. 
We clarify the concept of ergodic system giving a hierarchy of chaos. 



Let us consider a classical relativistic scalar field theory with action 

S[<P]= fd'xL{<P,d,<P), (1) 



where L is the Lagrangian density of the system, d^ = (^, V) is the covariant 
derivative and = 0(a;) is a real scalar field with x^ = (t, x) the space-time 
position^. It is well known that by imposing the Hamilton's Least Action 
Principle 

6S[<j>] = , (2) 

we obtain the Euler-Lagrange equation of motion of the system 

The homogenous space approximation means that we can neglect the spatial 
dependence of the field, thus we can perform the following substitution: 

0(t,x)-^0(t), (4) 

and the resulting equation of motion is given by 

dL d dL 

'dd>~Jt7hb~ ■ ^ ^ 



By introducing the momentum 






and the Hamiltonian 

H{xA) = h-L{<PA), (7) 

the second order equation of motion can be written as a system of two first 
order Hamilton's equations 

z = f (z) (8) 



where z = (0, x) is a point in a two dimensional phase space and f = (/i, /2) 

is given by: 

r)]-f f)M 

/i(0,x) = ^, /2(0,X) = -^. (9) 

This is a general result: any homogenous field theory can be written as 
a system of A^ first order differential equations, i.e. a dynamical system. In 
the next Sections we shall consider non-conservative and non-Abelian field 
theories. 

II-A Dynamical System Theory 

A dynamical system is defined by N first order differential equations 

z(t)=f(z(t),t), (10) 

where the variables z = {zi, ...,zn) are in the phase space Q (the euclidean 
space i?^, unless otherwise specified). These equations describe the time 
evolution of the variables and the system they represent^"^^. 

A solution of the dynamical system is a vector function z(zo,t), that 
satisfies (10) and the initial condition 

z(zo,0) = Zo. (11) 

Usually one writes simply z(t) without the initial condition dependence. 

The time evolution of z G ^2 is obtained with the one parameter group of 
diffeomorfism g^: Q ^ Q, such that 

|(^*z)|,=o = f(z,0). (12) 



The group g* is called phase flux and the solution is called orbit. The system 
is called Hamiltonian, if the dimension of Vt is even and there exists a function 
H{z,t) given by 

{{z{t),t) = JVH{z,t), (13) 



where: 



/ 

^H_, „l (14) 



is the symplectic matrix and H{z,t) is the Hamiltonian function. 

On the phase space Q one usually defines a probability measure fi : Q ^ 
Q, such that /i(fi) = 1. If we choose a subspace A of Q, the system is measure 
preserving if 

fiig'A) = fi{A) . (15) 

We observe that for measure preserving dynamical systems one gets div(f ) = 
0. It is well known that Hamiltonian systems preserve their measure: the Li- 
ouville measure. Dynamical systems which do not preserve their measure are 
called dissipative, and usually have a measure contraction in time evolution. 

The dynamics of a system is called regular if the orbits are stable to in- 
finitesimal variations of initial conditions. It is called chaotic if the orbits are 
unstable to infinitesimal variations of initial conditions. Useful quantities to 
calculate this behaviour are the Lyapunov exponents, which give the stability 
of a single orbit. 

A vector of the tangent space TQ^ to the phase space Q in the position 
z is given by 

^(z) = hmA^ ^, 16 

s~*0 s 



where r(0) = z and r(s) G Q. The tangent space vectors are the velocity 
vectors of the curves on M; there are obviously N independent vectors. 
Now we can define the Lyapunov exponent 

A(z)=lim^ln|cu(t)|, (17) 

where uj{t) is a tangent vector to z(t) with the condition that |cj(0)| = 1. 

It can be demonstrated that the limit given by the previous equation 
exists for a compact phase space, and that it is metric independent. Fixing 
an orbit in the N dimensional phase space, there are N distinct exponents 
Ai,...,A7v, called first order Lyapunov exponents. If the orbit has positive 
Lyapunov exponents, it is chaotic. 

To characterize globally the chaoticity of a system, we can introduce the 
Kolmogorov-Sinai entropy, which is given by 

hKsif^) = I M^) E ^*(z) ' (18) 

^ A,>0 

with A subspace of Q and Aj Lyapunov exponents. The Kolmogorov-Sinai 

entropy is a very useful tool for showing chaotic behaviour in the region A. 

A system is called ergodic if the time average is equal to phase space 

average 

hm i /* dtfiMt)) = f rf/i(z)/(z) . (19) 

i^oo t Jo Jn 

Incidentally, as is well known, Boltzmann started from the "ergodic hypoth- 
esis" to obtain statistical mechanics of equilibrium. But ergodicity is not 
sufficient to reach an equilibrium state: one must consider mixing systems. 

In a mixing system, every finite element of the phase space occupies for 
t ^ oo the entire phase space Q; more precisely: VA, B dVt with /i(A) and 
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/i(5) ^ 0, 



H"^^-^^^^ <-) 



To have quantitative information of orbit separations, we must introduce 
K-systems (Kolmogorov), which are mixing systems with a positive metric 
entropy, i.e. hxs > 0. Such systems are typical chaotic systems. 

Among the K-systems, the most unpredictable ones are the B-systems 
(Bernoulli), which have the Kolmogorov-Sinai entropy equal to the entropy 
of every partition, i.e. hxs = h{Ai{0),^), VAj(O). 

II-B Hamiltonian dynamics 

Let us consider an Hamiltonian system with n degrees of freedom described 
by the Hamiltonian function H{z), where z = (gi, ...g„,pi, ..., j9„) so that 
the phase space is A^ = 2n dimensional. The Hamiltonian system is called 
integrable if there are A^ functions Fi = Fi{z) defined on Q in involution: 

" dF- dF- dF- dF- 
[F^,F,U = T.p7f^-jf^P = 0, V.,j (21) 

^^ dqk dpk dqk dpk 

and linearly independent. [ , ]pb are the Poisson brackets. 
For conservative systems we have Fi = H{z) and also 

dF- 

^ = [H,F,]pB = 0. (22) 

Because there are n constants of motion, every orbit can explore only the 
n dimensional manifold Qf = {z : Fi{z) = fi,i = l,...,n}. If Qf is 
compact and connected, it is equivalent to a n dimensional torus T" = 
{{Qi, ■■■,Qn) mod 27r}. There are n irreducible and independent circuits 
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7j on Qf and there exists a canonical transformation (p, q) -^ (P,Q), gen- 
erated by the function S'(q, P), such that 

''• = *"''''■'' (23) 

The Pi are called action variables and the Qi are called angle variables. The 
moments p and coordinates q are periodic functions of Q with period 2tt. 
The Hamiltonian depends only on action variables, i.e. H = H(P). 

Adding a small perturbation V^(P, Q) to an integrable Hamiltonian Hq(P), 
the total Hamiltonian can be written: 

HiP,Q) = HoiP)+gViP,Q), (24) 

and, generically, the integrability is destroyed. As a consequence, parts 
of phase space become filled with chaotic orbits, while in other parts the 
toroidal surfaces of the integrable system are deformed but not destroyed; 
thus we have a quasi-integrable system. By growing g, chaotic motion de- 
velops near the regions of phase space where all the frequencies on the torus 
^i = Qp. are commensurate. Conversely, tori of the integrable system, on 
which the Ui are incommensurate, are deformed, but not destroyed immedi- 
ately (Kolmogorov-Arnold-Moser (KAM) theorem)^'^^. As g increases, the 
phase space generically develops a highly complex structure, with islands of 
regular motion (filled with quasi-periodic orbits) interspersed in regions of 
chaotic motion, but containing in turn more regions of chaos. As g grows 
further, the fraction of phase space filled with chaotic orbits grows until it 
reaches unity as the last KAM surface is destroyed. Then the motion is com- 
pletely chaotic everywhere, except possibly for isolated periodic orbits^'^^. 



It is very useful to plot a 2n— 1 surface of section V G il, called Poincare 
section. For an integrable system with two degrees of freedom, the gi = 
Poincare section of a rational (resonant) torus is a finite number of points 
along a closed curve, while the section of an irrational (non resonant) torus 
is a continuous closed curve. Adding a perturbation, the section presents 
closed curves (KAM tori), whose points are stable (elliptic), and also curves 
formed by substructures, residua of resonant tori, whose points are unstable 
(hyperbolic). As the perturbation parameter increases, the closed curves are 
distorted and reduced in number. 

II-C Quantum Chaos 

We use the term quantum chaotic system in the precise and restricted sense 
of a quantum system whose classical analogue is chaotic. In particular we 
concentrate on energy levels of quantum systems (see, for example, Ref. 18 
and 19). 

Let us consider a classical regular Hamiltonian system. The short-range 
properties of the corresponding quantal spectrum tend to resemble those of 
a spectrum of randomly distributed numbers. This is because regular clas- 
sical motion is associated with integrability or separability of the classical 
equations of motion. In quantum mechanics the separability corresponds to 
a number of independent conserved quantities (such as angular momentum), 
and each energy level can be characterised by the associated quantum num- 
bers. Superimposing the terms arising from the various quantum numbers, 
a spectrum is generated like that of random numbers, at least over short 



intervals. In particular, the distribution P{s) of nearest-neighbour spacings 
Si = (ej+i — ei)/d, where d is the mean level spacing, is expected to follow 
the Poisson limit, i.e. P{s) = exp (— s). 

Instead, when the classical dynamics of a physical system is chaotic, the 
system cannot be integrable and there must be fewer constants of motion 
than degrees of freedom. Quantum mechanically this means that once all 
good quantum numbers due to obvious symmetries etc. are accounted for, 
the energy levels cannot simply be labelled by quantum numbers associ- 
ated with certain constants of motion. The short-range properties of the 
energy spectrum then tend to resemble those of eigenvalue spectra of ma- 
trices with randomly chosen elements and one gets a result very close to 
P{s) = (7r/2)sexp (— ^s^), which is the so-called Wigner distribution. 

The distribution P{s) is the best spectral statistics to analyze shorter 
series of energy levels and the intermediate regions between order and chaos. 
This distribution can be compared to the Brody distribution 

P{s,u) = a{uj + l)s'^exp{-as'^+^) , (25) 

with 

a = (r[^])-^^ . (26) 

The Brody distribution interpolates between the Poisson distribution [uj = 0) 
of integrable systems and the Wigner distribution (u; = 1) of chaotic ones, 
and thus the parameter u can be used as a simple quantitative measure of 
the degree of chaoticity. 
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Ill A model for inflationary cosmology 

In this Section we study the stabihty of a scalar inflaton field^^'^^ and analyze 
its bifurcation properties in the framework of the dynamical system theory. 

It is generally believed that the universe, at a very early stage after the 
big bang, exhibited a short period of exponential expansion, the so-called 
inflationary phase. In fact the assumption of an inflationary universe solves 
three major cosmological problems: the flatness problem, the homogeneity 
problem, and the formation of structure problem^. 

The Friedmann-Robertson- Walker metric of a homogeneous and isotropic 
expanding universe is given by 

flv 

ds" = de - a\t) [ ^ + r\d9^ + sin^ 9d^% (27) 

where k = 1, —1, or for a closed, open, or flat universe, and a{t) is the scale 
factor of the universe. 

The evolution of the scale factor a{t) is given by the Einstein equations 

a = — —G{p + 3p)a, 



,d,2 k Stt 
a a'^ 6 



^-2= —Gp. (28) 



where p is the energy density of matter in the universe, and p its pressure. 
The gravitational constant G = M~^ (with h = c = 1), where Mp = 1.2 ■ 10"^^ 
GeV is the Plank mass, and H^ = a/ a is the Hubble "constant", which in 
general is a function of time. 

The inflationary models postulate the existence of a scalar fleld 0, the 
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so-called inflation field, with Lagrangian 

L = i9^09>-\/(0) (29) 

where the potential V{(f)) depends on the type of inflation model considered. 
The scalar field, if minimally coupled to gravity, satisfies the equation 



,a. ■ 1 „9 , dV 

- 0- — 
a a^ 



□0 = + 3(^)0 _ J_v20 = -^, (30) 



where □ is the covariant d'Alembertian operator. The Hubble "constant" 
Hu is related to the energy density of the field by 






Hl + ^ = & + ^ = ^[^ + ^^ + Vm. (31) 



In a fiat universe k = and, if the infiaton field is sufficiently uniform (i.e. 
0^, (V0)^ << V {(/))), we obtain an homogenous field theory in one dimension 

dV 
(P + 3H^{(P)(P + — = 0, (32) 

where the Hubble " constant" Hu is an explicit function of 0: 

Hi = ^V^(0). (33) 

III-A Local instability for the inflationary self— energy 

The second order equation of motion of our cosmological model can be writ- 
ten as a system of two first order differential equations 

z = f (z) (34) 
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where z = (0, x) is a point in the two dimensional phase space and f = (/i, /2) 
is given by 

/i(0,x)=X, /2(0,x) = -3//„(0)x-^^. (35) 

The system is non-conservative because the function 

is not identically zero. The fixed points of the system are those for which 
/i(0, X) = and /2(0, x) = 0, i.e 

The deviation 6z{t) = z(t) — z(t) from the two initially neighboring trajec- 
tories X and X in the phase space satisfies the linearized equations of motion 

j6z{t) = T{t)6z (38) 

where r(t) is the stability matrix 

/ 1 \ 

r(t) = , . (39) 

At least if an eigenvalue of r(t) is real the separation of the trajectories 
grows exponentially and the motion is unstable. Imaginary eigenvalues cor- 
respond to stable motion. In the limit of time that goes to infinity, from the 
eigenvalues of the stability matrix we can obtain the Lyapunov exponents. 
For two-dimensional dynamical system the Lyapunov exponents can not be 
positive^ and so the system is not chaotic, i.e. there is not global instability. 
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However, we can be assured that the universe is crowded with many inter- 
acting fields of which the infiaton is but one. The nonhnear nature of these 
interactions can result in a complex chaotic evolution of the universe and 
the local instability of the infiaton field is a precursor phenomenon of chaotic 
motion. 

The eigenvalues of the stability matrix are given by 



.,, = -|i/„W±i^9//2W-4^-12x^^ (40) 

The pair of eigenvalues become real and there is exponential separation of 
neighboring trajectories, i.e. unstable motion, if 

Particularly when x = 0, e.g. the fixed points, we obtain local instability 
when 

i.e. for negative curvature of the potential energy. The fixed points are stable 
if they are point of local minimum of V{(f)) and unstable if are points of local 
maximum. 

The potential V{(f)) depends on the type of infiation model considered, 
and it is usually some kind of double-well potential. We choose a symmetric 
double-well potential 

V{<j>) = ^{<P'-vy, (43) 

where ±f are the values of the infiaton field in the vacuum, i.e. the points 
of minimal energy of the system. 
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We observe that the inflaton field value in the vacuum f is a bifurcation 
parameter. Bifurcation is used to indicate a qualitative change in the features 
of the system under the variation of one or more parameters on which the 
system depends. First of all we consider the case v = 0, i.e. V{(f)) = (A/4)0^. 
In this situation there is only one fixed point (0* = 0, x* = 0) which is a 
stable one being 

^ = 3A0^ > 0. (44) 

The fixed point (0* = 0,x* = 0) is a point attractor. 
Instead for v ^ there are three fixed points 

(0* = O,x* = O), (0*=t;,x* = O), (0* = -t;,x* = O), (45) 

and the condition for the instability becomes 

V V , ^ 

Obviously (0* = 0,x* = 0) is an unstable fixed point, and in particular a 
saddle point because the stability matrix has real and opposite eigenvalues. 
On the other hand (0* = ±v, x* = 0) are stable fixed points. 
There are four possible functions for the Hubble "constant" 

HM = ±l\<P'-v% (47) 

but also 

HM = ±li<P'-v'), (48) 



where 7 = J2-kG\/3 is the dissipation parameter. The choice of the Hubble 
function is crucial for the dynamical evolution of the system. 
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In certain non-conservative systems we could find closed trajectories or 
limit cycles toward which the neighboring trajectories spiral on both sides. 
It is sometime possible to know that no limit cycle exist and the Bendixson 
criterion^^, which establishes a condition for the non-existence of closed tra- 
jectories, is useful in some cases. Bendixson criterion is as follows: if div(f ) 
is not zero and does not change its sign within a domain D of the phase 
space, no closed trajectories can exist in that domain. In our case we have 
div(f) = —3Hu{(f)), and so the presence of periodic orbit is related to the sign 

If Hu{(t)) = 7|0^ — f ^1 we do not find periodic orbits and the inflaton field 
goes to one of its two stable fixed points, which are points attractors (see 
Figure 1). The vacuum is degenerate but if we choose an initial condition 
around the saddle point there is a dynamical symmetry breaking towards 
the positive v or negative —v value of the infiaton field in the vacuum. This 
symmetry breaking is unstable because neighbour initial conditions can go 
in different point attractors. 

Instead, if we choose -f^„(0) = 7(0^ — v"^) the numerical calculations of 
Figure 2 show that exists a limit cycle, the two stable fixed points are not 
point attractors, and the infiaton field oscillates forever. Obviously more 
large is v more large is the limit cycle. 

III-B A limit cycle in the cosmological model 

The equation of motion of the infiaton field with if„(0) = 7(0^ — f^) reads 

+ 37(0^ - v'^)^ + A0(0=^ - t;2) = . (49) 
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This equation can be written as 

^ [0 + 37 f^iu^ - v^)du] + A0(02 - t;2) = , (50) 

at Jo 

and if we put 

F{(f))=3 f'^{u^-v^)du = (f){(f)^-3v^) , G{(f)) = (f){(f)^ - v^) , (51) 
Jo 

and also uj = (p + '~fF{(f)), we obtain the system 

^ ' ^^' (52) 

cl; = -AG'(0) . 

For systems of this kind the Lienard theorem^^ states that there is an unique 

and stable limit cycle if the following conditions are satisfied: -F(0) is an 

odd function and -F(0) = only at = and = ±0;; -F(0) < for 

< < a, -F(0) > and is increasing for > a; G'(0) is an odd function 

and 0^(0) > for all > a. It is easy to check that the functions -F(0) 

and G{(j)) satisfy all the conditions of the Lienard theorem with a = v. The 

cubic force (7(0) tends to reduce any displacement for large |0|, whereas the 

damping -F(0) is negative at small |0| and positive at large |0|. Since small 

oscillations are pumped up and large oscillations are damped down, it is not 

surprising that the system tends to Seattle into a self-sustained oscillation of 

some intermediate amplitude. 

Let us consider a typical trajectory of the system. After the scaling 

■0 = \uj we obtain 

= A[^-^F(0)], ^^^^ 

^ = -G(0). 
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The cubic nullcline ip = {'~f/X)F{(f)) is the key to understand the motion. 
Suppose that A >> 1 and the initial condition is far from the cubic null- 
chne, then we have |0| ~ 0{\) » 1; hence the velocity is enormous in the 
horizontal direction and tiny in the vertical direction, so trajectories move 
practically horizontally. If the initial condition is above the nullcline then 
> 0, thus the trajectory moves sideways toward the nullcline. However, 
once the trajectory gets so close that ip ~ (A/7)F(0) then the trajectory 
crosses the nullcline vertically and moves slowing along the backside of the 
branch until it reaches the knee and can jump sideways again. The period 
T of the limit cycle is essentially the time required to travel along the two 
slow branches, since the time spent in the jumps is negligible for large A. By 
symmetry, the time spent on each branch is the same so we have 

T ~ 2 / "^ rft (54) 

JtA 

where A and B are the initial and final points on the positive slow branch. 
To derive an expression for dt we note that on the slow branches with a good 
approximation ip ~ {'~f/X)F{(f)) and thus 

Since dip/dt = —0(0^ — f^), we obtain dt ^ ^3^-^, on the slow branches. 
The slow positive branch begins at 0^ = 2'~fv/X and ends at 0^ = '-fv/X. 



Because 7 = J27tGX/3 we get T ^ 21n2J^. 



IV The homogenous SU(2) YMH system 

Now we study the suppression of classical chaos in the spatially homogenous 
SU(2) Yang-Mills-Higgs (YMH) system induced by the Higgs field^'^^'^^ We 
analyze also the energy fluctuation properties of the system, which give a 
clear quantum signature of the classical chaos-order transition of the system. 
The SU(2) YMH system describes the interaction between a scalar Higgs 
field and three non-Abelian Yang-Mills fields A" a = 1,2,3. The La- 
grangian density of the YMH system is given by 

L = i(D,0)+(D^0) - \/(0) - iF;,F''- , (56) 



where 



iD,<P) = d,<P-tgAlT'<l), (57) 



F;, = d,At - d,A; + ge'^'^AlAl , (58) 

with T^ = a^/2, b = 1,2, 3, generators of the SU(2) algebra, and where the 
potential of the scalar field (the Higgs field) is 

\/(0)=/i'|0P + A|0r- (59) 

We work in the (2+l)-dimensional Minkowski space (/i = 0, 1, 2) and choose 
spatially homogeneous Yang-Mills and the Higgs fields 

9,A^ = 9,0 = , ^ = 1, 2 (60) 

i.e. we consider the system in the region in which space fluctuations of flelds 
are negligible compared to their time fluctuations. 
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In the gauge Aq = and using the real triplet representation for the Higgs 
field we obtain 

L = 0^ + i(i? + AD - g-[\AlAl - \{A, . A,f+ 

+ {Al + A^)02 _ (^^ . 0)2 _ (^^ . 0)2] _ ^.(0) ^ (61) 

where = (0\ 0^, 03)^ ^^ _ (^i^ ^2^ ^3) ^^^ ^^ _ (^1^ ^2^ ^3)_ 

When /i^ > the potential V has a minimum at |0| = 0, but for /x^ < 
the minimum is at 

,, , r^ 



'""' V 4A 

which is the non zero Higgs vacuum. This vacuum is degenerate and after 
spontaneous symmetry breaking the physical vacuum can be chosen 0o = 
(0,0,f). If A\ = qi, A2 = q2 and the other components of the Yang-Mills 
fields are zero, in the Higgs vacuum the Hamiltonian of the system reads 

H=l{pl + pD + g^v\ql + ql) + ]^g\lql , (62) 

where pi = gi and p2 = ^2- Here w"^ = 2g'^v'^ is the mass term of the Yang- 
Mills fields. This YMH Hamiltonian is a toy model for classical non-linear 
dynamics, with the attractive feature that the model emerges from particle 
physics. 

IV-A From chaos to order in the YMH system 

The chaotic behaviour of the YMH system can be studied by using the Toda 
criterion of the Gaussian curvature of the potential energy^'^'^^. For our YMH 
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system the potential energy is given by 

V{qu q2) = 9^v\ql + ql) + ]^g\lql . (63) 

At low energy, the motion near the minimum of the potential, where the 
Gaussian curvature is positive, is periodic or quasi-periodic and is sepa- 
rated from the instability region by a line of zero curvature; if the energy 
is increased, the system will be, for some initial conditions, in a region of 
negative curvature, where the motion is chaotic. According to this scenario, 
the energy E^ of chaos-order transition is equal to the minimum value of the 
line of zero Gaussian curvature KG{qi,q2) on the potential-energy surface. 
For our potential the gaussian curvature vanishes at the points that satisfy 
the equation 

{2gh' + g'ql) {2g%' + g'ql) - AgW.ql = . (64) 

It is easy to show that the minimal energy on the zero-curvature line is given 
by: 

Ec = VminiKc = 0, gi) = 6g\^ , (65) 

and by inverting this equation we obtain Vc = {E/6g'^Y^^. Thus the curvature 
criterion suggest that there is a order-chaos transition by increasing the 
energy E of the system and a chaos-order transition by increasing the value 
V of the Higgs field in the vacuum. Thus, there is only one transition regulated 
by the unique parameter E/{g'^v'^). 

It is important to stress that the Toda criterion is not a fully reliable 
indicator of chaos^^. In fact, the local instability of the Toda Criterion does 
not necessarily imply the global one and the idea of an order-chaos transition 
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with a critical energy is not strictly correct. The Toda curvature criterion 
should therefore be combined with the Poincare sections, which are shown 
in Figure 3. The numerical results confirm the analytical predictions of the 
curvature criterion: with E = 10 and g = 1 we get the critical value of the 
onset of chaos Vc = {E/Qg'^y^'^ ~ 1.14. 

IV-B Spectral statistics of the YMH system 

In quantum mechanics the generalized coordinates of the YMH system satisfy 
the usual commutation rules [qk,Pi] = i^ku with /c, / = 1,2. Introducing the 
creation and destruction operators 

«fc = y-gfc + «y— Pfc , "'k = y 2'^k - id —Pk , (66) 

the quantum YMH Hamiltonian can be written 

^ = ^0 + ^g'v , (67) 

where 

Ho = u^aiCii + ajd2 + 1) , (68) 

V = ^iai + atyia2 + at)\ (69) 

with cj^ = 2(7^f ^ and [ak, O;^] = Ski, k,l = 1,2. 

We compute the energy levels of the YMH system with a numerical di- 
agonalization of the truncated matrix of the quantum YMH Hamiltonian in 
the basis of the harmonic oscillators (see also Ref. 16 and 17). If \n1n2 > 
is the basis of the occupation numbers of the two harmonic oscillators, the 



22 



matrix elements are 

< n^n^\Ho\nin2 >= u;{ni + na + 1)5„;„^5„;„, , (70) 

and 

< n'i4|K|nin2 >= —[^ni{ni - 1)6^>^^^_^ + \/(^i + l)(«i + 2)5„;„^+2+ 



+ (2ni+l)(5„,„Jx [^n^K - lK'^^,^.,+ ^J {n2 + l)(n2 + 2)5„.„^_,2+(2n2+l)54„J 

(71) 
The symmetry of the potential enables us to split the Hamiltonian matrix 
into 4 sub-matrices reducing the computer storage required. These sub- 
matrices are related to the parity of the two occupation numbers rii and 
722: even-even, odd-odd, even-odd, odd-even. The numerical energy levels 
depend on the dimension of the truncated matrix: we compute the numerical 
levels in double precision increasing the matrix dimension until the first 100 
levels converge within 8 digits (matrix dimension 1156 x 1156). 

We have seen previously that the most used quantity to study the local 
fluctuations of the energy levels is the distribution P{s) of nearest-neighbour 
spacings Sj of the energy levels. It is obtained by accumulating the number 
of spacings that lie within the bin (s, s + As) and then normalizing P{s) to 
unit. 

We use the first 100 energy levels of the 4 sub-matrices to calculate 
the P{s) distribution. In order to remove the secular variation of the level 
density as a function of the energy E, for each value of the coupling constant 
the corresponding spectrum is mapped into one which has a constant level 
density. 
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The Figure 4 shows the P{s) distribution of Brody for three different 
values of the Higgs vacuum v. The best fit Brody parameter u is obtained 
by using the nearest-neighbour spacings of the first 100 unfolded energy levels 
of the YMH system. There is Wigner-Poisson transition by increasing the 
value V of the Higgs field in the vacuum. Thus, by using the P(s) distribution, 
it is possible to give a quantitative measure of the degree of quantal chaoticity 
of the system. Our numerical calculations show clearly the quantum chaos- 
order transition and its correspondence to the classical one. 

V Conclusions 

We have seen that spatially homogeneous field theories can be studied as 
dynamical systems. After a brief review of the dynamical system theory, we 
have discussed two schematic models of field theory. 

First, we have considered the stability of a non-conservative scalar in- 
fiaton field. The value of the infiaton field in the vacuum is a bifurcation 
parameter which changes dramatically the phase space structure. The main 
point is that for some functional solutions of the Hubble "constant" the sys- 
tem goes to a limit cycle, i.e. to a periodic orbit. The infiaton field is not 
chaotic but its local instability can give rise to a complex chaotic evolution of 
the universe due to its nonlinear interactions with other fields. In the future 
it will be very interesting to study these effects which can perhaps lead to 
some observable implications like a fractal pattern in the spectrum of density 
fiuctuations. 

We have then analyzed the non-Abelian SU(2) Yang-Mills-Higgs system. 
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We have given an analytical estimation (confirmed by numerical results of 
Poincare sections) of the classical chaos-order transition as a function of 
the Higgs vacuum, the Yang-Mills coupling constant and the energy of the 
system. A quantum signature of a chaos-order transition has been obtained 
by using the distribution P{s) of nearest-neighbour spacings. The Wigner- 
Poisson transition of the P{s) distribution follows very well the classical 
results of the Poincare sections. 

To conclude, we observe that there are yet many open problems about 
chaos in field theory. We make a list of some of them: i) spatial chaos and 
space-time chaos; ii) classical and quantum chaos in more realistic systems, 
for example in QCD (some results can be found in Ref. 7 and 8); iii) connec- 
tion between chaos and critical phenomena (finite temperature field theory). 
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Figure Captions 

Figure 1: The Hubble function wstime (top) and tlie pliase space trajectory 

of tlie infiaton field (bottom); for if„(0) = 7|0^ — f^| with 7 = 1/2, A = 3 

and V = 1. Initial conditions: = and = 1/2. 

Figure 2: The Hubble function fstime (top) and the phase space trajectory 

of the infiaton field (bottom); for Hu{(f)) = 7(0^ — v"^) with 7 = 1/2, A = 3 

and V = 1. Initial conditions: = —1/2 and = 0. 

Figure 3: The Poincare sections of the YMH system. From the top: v = 1, 

V = 1.1 and V = 1.2. Energy E = 10 and interaction g = 1. 

Figure 4: P{s) distribution of Brody of the YMH system. First 100 energy 
levels and g = 1. The best fit Brody parameter is given by: ui = 0.92 for 

V = 1.0, 00 = 0.34 for V = 1.1 and uj = 0.01 for v = 1.2. 
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